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Abstract 

How to macroscopically control the flow of heat at will is up to now a challenge, which, however, 
is very important for human life since heat flow is a ubiquitous phenomenon in nature. Inspired 
by intelligent electronic components or intelligent materials, here we demonstrate, analytically and 
numerically, a unique class of intelligent bifunctional thermal metamaterials called thermal cloak- 
concentrators, which can automatically change from a cloak (concentrator) to a concentrator (cloak) 
when the applied temperature held decreases (increases). For future experimental realization, the 
behavior is also confirmed by assembling homogeneous isotropic materials according to the effective 
medium theory. The underlying mechanism originates from the effect of nonlinearity in thermal 
conduction. This work not only makes it possible to achieve a switchable Seebeck effect, but also 
offers guidance both for macroscopic manipulation of heat flow at will and for the design of similar 
intelligent multifunctional nretanraterials in optics, electromagnetics, acoustics, or elastodynamics. 
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I. INTRODUCTION 


The flow of heat is a common phenomenon in nature. However, how to macroscopically 
manipulate the heat flow at will remains poorly explored in the literature, which, however, is 
very important for human life or military uses, say, efficient refrigerators, solar cells, energy¬ 
saving buildings, or infrared camouflage. It is known that our world has been changing by 
intelligent electronic components or intelligent materials. Here the intelligent means that the 
component or material has properties that can be tuned in a controllable manner by external 
stimuli. For example, electronic diodes belong to intelligent electronic components, which 
have much lower resistance to electric current along one direction than the other; shape 
memory alloys are a kind of intelligent materials, which can recover their original shapes 
by stimuli (say, heating [1] or adding an external magnetic field ^|) after big deformations. 
Manipulating thermal conduction equations to achieve desired temperature distributions or 
heat flux patterns in thermal metamaterials is promising for thermal extensions of intelligent 
electronic components or intelligent materials. As a result, this manipulation might produce 
a unique class of intelligent thermal metamaterials, thus yielding a novel (intelligent) way 
for macroscopically controlling the flow of heat at will. For this purpose, let’s trace back to 
the birth of transformation optics theory 3, 4], 

In 2006, Leonhardt and Pendry et al. independently established the transformation optics 
theory, which leads to a new concept of invisibility cloaking [ 3 , 4], Such a significant progress 
soon enlightened a lot of scientists in different fields since it offers a powerful method for de¬ 
signing metamaterials with new properties like optical illusions js-7]. I 11 essence, this theory 
demands the form invariance of domain equations under linear coordinate transformations. 
This characteristic serves as the foundation of transformation optics theory. Inspired by it, 
some researchers extended the transformation mapping theory to the domain of thermal con¬ 
duction, and then proposed a type of thermal metamaterials, which are called thermal cloaks, 
for steady state heat flow (i.e., the temperature, T, is independent of time) by using inhomo¬ 
geneous anisotropic materials 8]. The original intention of designing such type of thermal 
metamaterials is to hide an object inside the cloak from the detection through measuring the 
external temperature distribution, and thus these metamaterials can thermally protect the 
central region containing the object. So far, the steady-state thermal cloak fs and its theo¬ 


retical extensions 


10] have been experimentally realized or developed llMl5l|. Meanwhile, 
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researchers also designed a lot of thermal metamaterials wit 


beyond cloaking 
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24], such as concentrators 
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lnovel thermal characteristics 
16]. The concentrator helps to 


focus heat flux in a particular region, and thus yields a higher temperature gradient inside 
the region, which can be used to improve the efficiency for thermal-to-electrical conversion 
(the Seebeck effect [25, 26]). 

Unfortunately, the existing thermal metamaterials published in the literature QQ don’t 
have the essential feature of intelligence, namely, they lack the ability to automatically 
sense and respond to the change of environmental temperature. The reason is that within 
the framework of transformation thermotics (or thermodynamics as also named by some 
other famous scholars) 


0, Q, 


the property of thermal metamaterials can influence the 
temperature held, but not vice versa due to temperature-independent thermal conductivities 
under consideration. Evidently, it becomes necessary to develop the theory of nonlinear 
transformation thermotics 24], which takes into account the fact that many materials have 


a temperature-dependent thermal conductivity. This theory suggests that new function can 
be introduced into thermal metamaterials by using a temperature-dependent transformation. 
In this work, we shall show that thermal metamaterials can be made intelligent indeed based 
on the above idea by presenting a new thermal metamaterial — intelligent thermal cloak- 
concentrator, which can automatically change from a cloak (concentrator) to a concentrator 
(cloak) when the applied temperature held decreases (increases). 


II. METHOD 


The traditional thermal cloak helps to prevent heat from traveling inside the central 
region without disturbing the temperature distribution outside, and thus thermally hides 
an object inside it. As schematically drawn in Fig. 1(a), a two-dimensional thermal cloak 
is designed according to a compression transformation which maps a circle of radius R 2 to 
a ring with interior radius Ri and exterior radius 1?2- In a polar coordinate system, the 
transformation between original coordinates (r, 9) and transformed ones (r', 6') is written as 




R-2 — R\ 

Ro. 


+ R 1 ,9' = 9. 


( 1 ) 
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According to this transformation, the thermal conductivity k of the cloak is calculated from 
the equation of traditional transformation thermotics, 


J no R 


K det (J)’ 

where J is the Jacobian matrix of the transformation written in Cartesian coordinate system, 
is the thermal conductivity of the background, and det(J) is the determinant of J. 

The traditional thermal concentrator guides more heat flux into the central region and 
increases the temperature gradient inside. As Fig. 1(b) shows, the corresponding transfor¬ 
mation of a two-dimensional concentrator maps a ring with interior radius R 3 and exterior 
radius f? 2 to a ring with interior radius R\ and exterior radius Ro. The formulas of this 
transformation in a polar coordinate system are 

-,0' = 9 


/ R 2 — R\ , n R\ — R 3 




+ R2 


( 3 ) 


R2 — R3 R2 — R3 

The thermal conductivity of the concentrator can also be calculated by substituting the 
Jacobian of Eq. ((3]) into Eq. ([2]). 

A comparison between Eq. (JT]) and Eq. (J3J) reveals an interesting fact: if we replace the R 3 
in Eq. (J3J) with 0, we get exactly the transformation in Eq. (IT]) . The reason behind it is that 
the original configuration of the transformation is a circle, and can be regarded as a ring with 
interior radius 0. This resemblance between the principle of the cloak and the concentrator 
inspires us to combine these two thermal metamaterials together and design a new kind of 
intelligent materials. This device will automatically switches from a concentrator to a cloak 
or vice versa as environmental temperature changes. However, the traditional transformation 
thermotics does not serve as a proper tool to design it because the temperature won’t affect 
the behavior of the device. Thus an extension of the traditional theory is needed. 

In our previous work 24], we’ve shown that the transformation thermotics can be gener¬ 
alized by introducing the nonlinear effect of temperature-dependent thermal conductivities. 
This theory enables us to make use of a certain temperature-dependent transformation to 
achieve novel functions such as switchable thermal cloaking and macro thermal diodes based 
on the cloaking 24]. The foundation of this approach is the equivalence between such trans¬ 
formation and a normal transformation applied on nonlinear materials (whose conductivity 
depends on temperature). 

Since we can change a concentrator into a cloak by simply replacing R 3 in Eq. (J3J) with 
0, the cloak-concentrator can be designed on the basis of a function R*(T ) which switches 
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between R% and 0 at different temperature T. Then the T-dependent transformation is 




R 2 - R\ + r ^ R\ ~ R*(T) q, _q 


( 4 ) 


R 2 -R*(T) R 2 — R*(T) 

The thermal conductivity k(T) of our device is calculated from a relation similar to Eq. (J2]), 

K(T) = det Jim ’ (5) 

where J(T) is the Jacobian matrix corresponding to the transformation in Eq. ([4j) written 
in the Cartesian coordinate system. The result expressed in the polar coordinate system is 




1 + 


R 2 (R*(T) — Ri 
r'(R 2 — R*(T)) J 


,Ko(T) = Kq 


1 + 


RviR'W-Rd 

r'(R 2 -R*(T )) 


-1 


( 6 ) 


where K r (T) and Kg(T) are the radial and azimuthal components of k(T), respectively. 

The specific term of function R*(T) is determined according to the desired behavior of 
the cloak-concentrator. In this study, we want the device to behave like a concentrator when 
the environmental temperature is lower than a critical temperature, T c , such that the central 
region inside it receives more heat flux from the heat source. And when the environmental 
temperature increases to be higher than T c , we want the device to behave like a cloak in 
order to avoid heat flux from coming inside the central region. In this case, an object located 
in the central region is able to experience a relatively more stable thermal environment. To 
achieve the above function, we resort to the following form of R*(T), 


R*(T) = R 3 / [1 + e^ (r - Tc) ] , 


(7) 


where /3 is a scaling factor which determines the sensitivity of the device to temperature. 


III. RESULTS 

In order to show the effect of the bifunctional device, we perform finite element simula¬ 
tions based on the commercial software COMSOL Multiphysics to see if the metamaterials 
work when the applied temperature fields change. Without loss of generality, the units of 
length, temperature and thermal conductivity are defined as L x /8 (L x is the width of the 
simulation box), the initial temperature and the conductivity of background, respectively. 
Accordingly, the scales of the length, temperature and conductivity are nondimensionalized 
throughout this work. An annulus with interior radius R,\ — 1 and exterior radius R 2 = 3 
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is set in a box with size 8 x 7 as shown in Fig. 2. The thermal conductivities of the back¬ 
ground and the central region are set as Kq = 1. The thermal conductivity of the device is 
calculated according to Eq. ((0]) for R 3 = 2. Heat diffuses from the left boundary with high 
temperature T B to the right boundary with low temperature T L . Meanwhile, the upper and 
lower boundaries of the simulation box are thermally isolated. 

The simulation results of the switchable device are shown in Fig. 2(a,b) where the critical 
temperature T c is set to be 1.5. In Fig. 2(a) for low temperature (T = 1.2 ~ 1.4), we observe 
that the concentrator is focusing the heat flow into the central region and the temperature 
gradient is raised. On the other hand, when the environmental temperature is high enough 
(T = 1.6 ~ 1.8) as in Fig. 3(b), the cloaking effect is “turned on”. The heat flux in the 
cloak are guided to detour around the core (namely, the central region) and the temperature 
inside this core is turned to be a constant. Due to the antisymmetry, a switchable device 
between a cloak in low temperature and a concentrator in high temperature can also be 
worked out on the same footing. 

Although the thermal conductivity of such a device is inhomogeneous and anisotropic 
[Fig. 2(a,b)], for future experimental realization we may resort to the effective medium the- 

m 


ory (EMT) 


28] which utilizes a laminated structure as an approximation of anisotropic 


materials. According to the EMT, one layer contains two alternating homogeneous isotropic 
sub-layers of thicknesses and d B with conductivities Ka(T) and k b (T). In a polar coor¬ 
dinate system, the effective parameters for a layer can be obtained from 


+ 


V 


,ko(T) = 


k a (T) + tjkb{ T) 

1 + r] 


( 8 ) 


K r (T) 1 + T)\I%a(,T) K B (T ) 

where r/ = d B /dA- In this work, we set 77 = 1 for simulations. 

The method based on the EMT has been successful in the experimental demonstration 
of thermal cloaks jlT, 12], For the case of concentrators, the transformation is based on 
extensions; therefore, the array orientation of the alternating sub-layers should be rotated 
with an angle of 7 t/ 2. Considering both the aforementioned transformations, as shown in 
Fig. 2(c,d), we divide the device into 400 grids (which contain 5 layers consisting of 10 
sub-layers in the radial direction and 20 layers consisting of 40 sub-layers in the azimuthal 
direction). 

Eq. (JED is considered to be a reliable tool in realizing the theoretical design, but the 
resulting forms of ka(T) and k b (T) obtained by solving Eq. (JED are too complicated. So 
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in this study, we propose that instead of getting the thermal conductivities directly from 
Eq. ()8j). we solve two limiting cases of concentrators and cloaks separately, 



1 


1 ( 1 


.-Cl I nr) Lc Cl 

, Cl _ ^ + r l K B 

— -i i 


( 9 ) 


In the above equations, the superscript, Co (or Cl), denotes the case of traditional concen¬ 
trators (or cloaks) without temperature dependency. Then, according to the solutions of 
Eq. ([9]), we find the desired relations of ka(T) and kb(T) with respect to temperature, 


ka(b){T) — «5(b) + { k a°b) ~ k a(b )) [l + e/3(T Tc) ] 


( 10 ) 


Now we are in a position to assemble the facility. Evidently, Fig. 2(c,d) displays the same 
switching effect between concentrators and cloaks as presented in Fig. 2 (a, b). 

Furthermore, for clarity, in Fig. 3, we choose the temperature gradient inside the device 
to characterize the state of the intelligent thermal metamaterial. Since q = kVT, where 
q is the heat flux and n is the thermal conductivity of the central region (k = 1), the 
temperature gradient represents the amount of heat flux in ^-direction. A high value of 
gradient in the central region indicates the device is working as a concentrator, while the 
near zero value denotes that the cloaking effect is activated. For different temperature fields 
applied, the variation of the temperature gradients are displayed in Fig. 3. The results show 
that the new thermal metamaterial can make a quick response to the change of environmental 
temperature, and the transition process is sharp around the critical point T c ( T c = 1.5). The 
simulation results based on the EMT are also provided and the behaviors of the metamaterial 
are in good agreement with the theory. The slight amount of temperature gradient at the 
cloak side is due to a relatively few number of layers (5 layers as plotted in Fig. 3) in the 
radial direction. This non-zero temperature gradient can be further reduced to approach zero 
with more alternative rings of thinner thickness. As an evidence, the concentrator effects 
for the nonlinear transformation theory and EMT are very close to each other because we 
used 20 layers of materials in the azimuthal direction (Fig. 3). 
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IV. DISCUSSION AND CONCLUSION 


By taking into account the role of nonlinearity in thermal conduction, we have designed a 
unique class of intelligent bifunctional thermal metamaterials — thermal cloak-concentrators 
with intelligence. The cloak-concentrator is able to function as two traditional metamaterials 
(namely, a cloak and a concentrator) at different environmental temperatures and quickly 
switches from one to the other as temperature changes. With this new device, the heat flux 
amount of a certain region can be automatically adjusted according to the environmental 
temperature. For the convenience of future experimental demonstration, we have performed 
simulations based on the design of assembling homogeneous isotropic materials according to 
the EMT (effective medium theory); see Fig. 2(c,d). The thermal conductivities of theses 
materials vary with temperature as determined by Eq. m, indicating that the value of 
conductivities should increase or decrease rapidly with temperature. Fortunately, researchers 


have found several kinds of materials whose thermal conductivities have this behavior 


M 


30], which provide a viable approach to realize the intelligent thermal cloak-concentrator 


proposed in this work. 

Owing to the compression transformation of the original coordinate system, the perfect 
cloaking effect provides a dead core inside the device, where the temperature gradient are re¬ 
duced to zero. On the contrary, the concentrator can raise the gradient in the central region 
due to the extension transformation of the shell. Since for a given thermal conductivity the 
temperature gradient determines the amount of heat flux, the thermal cloak can significantly 
reduce the heat flux inside it, but the thermal concentrator can raise the heat flux amount 
inside to an abnormally high value. Naturally, if one device can switch between a cloak and 
a concentrator, the amount of its inner heat flux can vary in a wide range; if the switching 
effect can be triggered by changing the environmental temperature, the desired intelligence 
is achieved. Such a bifunctional thermal metamaterial can automatically cause its central 
region to heat (or cool) in case of a low (or high) environmental temperature. In particular, 


it can be used to switch on or off the Seebeck effect 


25 


26] when the applied temperature 


varies, thus achieving a so-called switchable Seebeck effect. Besides, it provides guidance 


for macroscopically manipulating heat flow at will in many thermal management 


which are related to heat preservation, dissipation, camouflage 


20] or illusions 


pro 


22 


Diems, 
23J in 


various kinds of areas ranging from energy-saving machines/buildings to military uses. In 
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addition, due to the similarity between the governing equations, our consideration adopted 
in this work can be extended to obtain the counterparts of such intelligent multifunctional 
thermal metamaterials in other disciplines like optics, electromagnetics, acoustics, and clas- 
todynamics. 
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Figure Captions 

Fig. 1. Schematic graph showing (a) a thermal cloak and (b) a thermal concentrator. 
Heat flow is represented by red lines with arrows. The interior radius and exterior radius of 
the device is R 1 and i? 2; respectively. The area with radius Ri denotes the central region or 
core. 

Fig. 2. Finite-element simulations of the cloak-concentrator. (a,c) The device functions 
as a concentrator at temperature below 1.4 and (b,d) functions as a cloak at temperature 
higher than 1.6. The temperature distribution is denoted by the color, with the white 
isothermal lines indicated. The left and right boundaries are heat and cold sources with 
fixed temperatures. And the upper and lower boundaries are thermally insulated, (a) and 
(b) show the results of Eq. (EJ); (c) and (d) display the results of Eq. (ITU]) . Parameters: 
Ri = 1, i ?2 = 3, A 3 = 2, kq = 1, and T c = 1.5. 

Fig. 3. Relation between the environmental temperature (defined as the average temper¬ 
ature of the heat and cold sources on the left and right boundaries) and the temperature 
gradient of the central region (which corresponds to the area of r < R\ in Fig. 1) along x 
direction. As shown in the figure, the red circles and blue squares denote the results based 
on the nonlinear transformation theory [Eq. (J 6 ])] and EMT [Eq. (fTOl) ] respectively. The re¬ 
direction gradient outside the device keeps constant, and it is indicated by the horizontal 
dashed line; the critical temperature T c is represented by the vertical dashed line. 
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Temperature Gradient 
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